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Shear viscosity is evaluated within a model of the gluon plasma, which is based 
entirely on the stochastic nonperturbative fields. We consider two types of excitations 
of such fields, which are characterized by the thermal correlation lengths ~ (g^T)^^ 
and ~ (g'^T)^^, where g is the finite-temperature Yang-Mills coupling. Excitations 
of the first type correspond to the genuine nonperturbative stochastic Yang-Mills 
fields, while excitations of the second type mimic the known result for the shear 
■ viscosity of the perturbative Yang-Mills plasma. We show that the excitations of the 

^ first type produce only an O (g^^)- correction to this resuh. Furthermore, a possible 

(— 1 1 interference between excitations of these two types yields a somewhat larger, 0{g'^), 

correction to the leading perturbative Yang-Mills result. 

Our analysis is based on the Fourier transformed Euclidean Kubo formula, which 
represents an integral equation for the shear spectral density. This equation is solved 
by seeking the spectral density in the form of the Lorentzian Ansatze, whose widths 
are defined by the two thermal correlation lengths and by their mean value, which cor- 
responds to the said interference between the two types of excitations. Thus, within 
one and the same formalism, we reproduce the known result for the shear viscos- 
ity of the perturbative Yang-Mills plasma, and account for possible nonperturbative 
corrections to it. 



I. INTRODUCTION 



Over the last ten years, it has been widely discussed that the quark-gluon plasma pro- 
duced in the RHIC experiments can resemble an almost perfect quantum liquid, which is 
characterized by the^values of the shear- viscosity to the entropy-density ratio, r]/s, much 
smaller than unity [1]. Comparison with the empirical data for water, helium, and nitro- 
gen shows that their {rj/ s)-Tatios reach minima in the vicinity of the corresponding liquid- 



gas phase transitions Given different types of phase transitions and different types of 
molecules for the above-mentioned three substances, one can guess that such a temperature- 
behavior of rj/s is quite general. Using this observation, one can naturally assume that for 
the quark-gluon plasma the minimum of the {ri/s)-iaiio is reached in the vicinity of the 
deconfinement phase transition. This minimum can be set to the value of l/(47r), obtained 
within {M = 4) supersymmetric Yang-Mills theory, which was suggested as the lower bound 
for the (?7/s)-ratio jsl. With the increase of temperature, rj/s is expected to rise from this 
bound up to the values predicted by the perturbative-QCD calculations j4|. Thus, a problem 
can be posed as how to model these essential features in the temperature-behavior of the 
{ri/s)-T8itio. 

In the present Letter, we address this issue for the purely gluonic plasma, within a 
model based entirely on the stochastic nonperturbative fields. In particular, we manage 
to reproduce the said high-temperature behavior of rj/s in perturbative Yang-Mills theory 
by imposing the correlation length of these fields to be ~ {g^T)~^, where g is the finite- 
temperature Yang-Mills coupling. This correlation length is recognizable as the mean time 
needed for a parton undergoing Coulomb scatterings in the gluon plasma to deflect by an 
angle of the order of unity [5]. Of course, besides the ultrasoft momentum scale ~ g^T, 
stochastic nonperturbative fields possess just the soft scale ~ g^T, which defines the high- 
temperature behavior of such quantities as the spatial string tension [6] and the nonpertur- 
bative gluonic condensate Hence, the key ingredient of our model is the presence, in the 
deconfinement phase (T > Tc) of interest, of the two types of excitations of the nonpertur- 
bative fields. These excitations are characterized by the parametrically different correlation 
lengths, ~ {g'^T)~^ and ~ (g'^T)^^. Excitations of the first type describe an extrapolation 
of the genuine nonperturbative stochastic vacuum Yang-Mills fields to the deconfinement 
phase |7|. Rather, excitations of the second type are introduced with the purpose to mimic 
the known perturbative contribution Q| to the shear viscosity. 

The goal of the present study is therefore twofold: to quantify the relative contribution to 
rj/s, which is produced by the genuine nonperturbative fields [i.e. those with the correlation 
length ~ {g^T)^^ with respect to the known perturbative contribution, and to evaluate the 
contribution to r]/s produced by the perturbative-nonperturbative interference. These two 
issues will be addressed by obtaining the shear spectral density from an integral equation 
given by the Fourier transformed Kubo formula. That will be done by seeking the spectral 
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density as a superposition of the Lorentzian Ansatze. As a result, we find that the widths 
of these Lorentzians are given by the said momentum scales, ~ g'^T and ~ g'^T, as well 
as by their mean value (for the interference between the perturbative and nonperturbative 
interactions). 

In the next Section, we perform the corresponding analytic and numerical calculations. 
In Section III, we give a brief summary of the results obtained. 



II. CALCULATION OF THE (ry/s)-RATIO 

The spectral density p = p{u,T), defining the shear viscosity t] = rj(T) as 

dp 

1 

can be obtained from the following Euclidean Kubo formula jsl: 



7] = vr — 

auj 



(1) 



smh(a;/3/2) J ^ 

^ ' ' n=— oo 

Here, P = 1/T, n labels the winding mode, and Ut is the finite-temperature correlation 
function of the (1, 2)-component of the Yang-Mills energy-momentum tensor 0^i,: 

f/r(x,X4 + /3n) = (612(0, 0)0i2(x,X4 + /3n)>^, where = g^F^^F^^. (3) 

The Kubo formula represents an integral equation for p. To solve this equation, we first 
Fourier transform it. This method of solving the equation is inspired by the observation 
that 

cosh [g; (X4 - f )] _ ^ e-'--^ 

sinh(c./3/2) ■^^Z^^u^^ + ujr 

where Uk = 2'KTk is the k-th Matsubara frequency. One further notices that, for non- 
perturbative fields at issue, Ut exponentially falls off at a distance defined by the thermal 
correlation length of those fields. For this reason, we consider the maximally general expo- 
nential Ansatz for Ut-, which is provided by the MacDonald functions. Namely, we start 
with the following sum, which generalizes the one on the right-hand side of Eq. 
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Here, r(a) stands for the Gamma-function, a > 0, and m = m{T) is some mass parameter. 
The sum over Matsubara modes k can be transformed into a sum over winding modes n, 
yielding the following intermediate result: 



oo 

e 

n=— oo 



One can further multiply this expression by 1 = (47rA)^'^/^ J (Pxe~^, which yields 
Performing then the A-integration, one obtains 



-oo 

e 

n=— oo 



7r'r(a) 7 ^f-^ (mv/x2 + (^4 + f3nff- ' 

where Ki,{x) stands for the MacDonald function. 

Hence, we assume the correlation function ([Sj) of the following form: 

f/T(x, X4 + = Nia) Gt =— -, (5) 

(mA/x^ + (X4 + pnyy " 

where A^(a) is a numerical parameter, and = {{gF^y)'^)T is the finite-temperature non- 
perturbative gluonic condensate. Then the Fourier-transformed Kubo formula reads 

duj = 7r22"iV(a) T{a)Gl ^. (6) 

We use now for p a Lorentzian Ansatz with the width m: 

Cu 



where C = C{T) is the sought function of dimensionality (mass)^. Notice that, although 
the asymptotic freedom requires p oc a.t u ^ T (cf. e.g. Ref. 13|), it is the Lorentzian 
part of p which matters for r], since it defines the derivative of p at w = 0. As it then follows 
from Eq. ([T]), the shear viscosity is expressed in terms of C as rj = TtC/m?. Substituting 
Ansatz ([7]) into Eq. ([6]), we are left with the integral 



du 
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Setting a = 1/2, we arrive at the relation 



C = VS^N 



G"^ \uJk\ + m 



where from now on = N{l/2). Thus, we see that, for a = 1/2, an exponentially falling 
off function is compatible with the Lorentzian Ansatz ([7]) for A; = (that is, in the high- 
temperature dimensionally-reduced theory) and for \k\ ^ 1. With a given form (jl]) of the 
kernel in the integral equation ([2]), which is prescribed by the fluctuation-dissipation theorem, 
and with the use of the Lorentzian Ansatz, a better accuracy can hardly be achieved. Thus, 
we use the formula 

nc^Vs^N^, (8) 
which is supported by the observation that the correcting factor 

(9) 

is equal to 1 for T > T^,, where T^, is the temperature of dimensional reduction. By the end 
of our analysis, we will numerically evaluate maximum possible deviations of the correcting 
factor from 1, which can take place for the temperatures T G (Tc,T^,) at \k\ ~ 1. 

We proceed now to the calculation of the coefficient A^. To this end, we first notice that 
the (T = 0)-counterpart of the function ([5]) at a = 1/2 reads 

TT f \ ^3/2(mo|a;|) ^^3/2(2) 1 A, .A -z n n^ 

""'^'^ = (mo|x|)3/2 ' "^^^^ = V 2 U J ' ' ^ ^ 

G = {{gFp^), and the subscript "0" means "at zero temperature". Second, we express 
this correlation function in terms of the 2-point functions of -F^j^'s by using the so-called 
Gaussian-dominance hypothesis. This hypothesis, supported by the lattice simulations ^, 
states that the connected 4-point function of -F^^^'s can be neglected compared to the pairwise 
products of the 2-point functions. For the function Uq{x) = {g^F^fj_{0)F2^{0)F^^{x)F2^,{x)) 
at issue (cf. Eq. ([2])), the Gaussian-dominance hypothesis yields 

Uoix) ^ {9'F,^,{0)Fl{x)){g'F,^^{0)Fl{x)) + {9'F^,{0)Fl{x)){g'F,^^{0)Fl{x)), (11) 

where we have taken into account that {g'^ F^^{0) F2^{0)) = 0. The contribution of stochas- 
tic nonperturbative fields to the function (5'^-F^,y(0)-F^p(x)) can with a high accuracy be 
parametrized as 

m 

G S""^ 
{g^F;MnM)) = (^M^.p - 5,p5.A) ■ ■ Dix). (12) 
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In Eq. (fT2|) . the dimensionless function D{x) exponentially falls off at a distance called 
the vacuum correlation length. Substituting Eq. (fT2l) into Eq. f|TT]) . one readily obtains 
Uo — 72(jv2-i) -^^- Comparing this expression with Eq. (fTOl) . and setting from now on Nc = 3, 
we have 

D(x) = 24 



' ^ Ks/2imQ\x\) 



1/2 

(13) 



(mo|a;|)^/2 

As mentioned above, we consider two types of excitations of the nonperturbative fields, 
which are characterized by the correlation lengths ~ (g'^T)^^ and ~ {g'^T)~^. We start 
with the excitations of the first type. These excitations exist foremost in the confinement 
phase (T < T^), where they yield the string tension [10] cr = ^ f (PxDiy.) corresponding 
to the static sources in the fundamental representation. On the other hand, the function 
D{x) in the confinement phase is conventionally parametrized by just an exponent jol, [l^ 
which, once compared to Eq. ( IT3|) . would be D{x) = e~™°l^l/^. Plugging both this expo- 
nential parametrization and parametrization (|T3|) into the said formula for a, we obtain the 
corresponding coefficient N: 

N = ^ , (14) 

[6j^°^dxxy^,/K^]' 

From the confinement phase, we immediately jump to the opposite limit of very high 
temperatures, T ^ T^, where the excitations of the second type are supposed to be mostly 
important. There, the shear viscosity has the form [4] 



_ 27.126 



At such temperatures [l], Gt ~ (fi-^T)^, so that Eq. (ED yields A^' ~ (^)°. [We use the 
notation A^' to make a distinction from Eq. (IHj).] In order for this coefficient A^' to be 
constant, one should have m ~ g'^T. Thus, the known high-temperature expression for the 
shear viscosity can indeed be reproduced within a model of nonperturbative stochastic fields 
with the correlation length ~ {g'^T)~^. 

To distinguish the two scales, ~ g'^T and ~ g'^T, we use from now on the notations M 
and m, respectively. Thus, 

M = mofr and m = g^M. (16) 
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Here, the continuous function Jt can be chosen in the following form: 

, 1 1/4 



h 



The coth-f actors in Eq. ([T 



COthf 
COth( 



at Tr < T < 



1/4 



(17) 



at T > T^. 



stem from the above-mentioned parametrization D{x) = 
which yields H ^ = /| = coth(^f^j at < T < T,. Furthermore, T, in Eq. ([HD 
is the temperature of dimensional reduction, and = q{T^). Lattice simulations Igl and 

nn 

analytic calculations p, HH suggest the value of in the range from Tc to 2Tc. Below, we 
will find from this range by using the best known lattice value for the shear-viscosity to 
the entropy-density ratio. 

We assume now the function Dix) aX T > in the form of a sum 



D{x) =DM{x) + Dm{x). 



(18) 



Here, Dm is given by Eq. (fT3l) with tjiq replaced by M, while is given by a similar formula: 



Drr,.{x] 



24 



■ , K-j,/2(m\x\) 
(m|a;|)3/2 



1/2 



. Since M ^ m at T ^ T^, one has at such temperatures 
D ^ Dm with an exponentially high accuracy. Therefore, the shear viscosity ([H]) goes at 
T ^ Tc as 



r/ ~ N' 



coth( ^ 



1/4 



Comparing this expression with the known one, Eq. (fT5i) . we get the coefficient A^': 

3 



N' 



1 27.126 



VSttS In 



2.765 



COthf 



1/4 



m5 

^2 



As it was anticipated above, this result is T-independent with the double logarithmic accu- 
racy (since g depends on T only logarithmically). 

We can now proceed towards the main result of the present paper — the full shear 
viscosity produced by the two types of excitations of the nonperturbative stochastic fields, 
which are described by the correlation function (ITSl) . 



Dix] 



24 



A^ 



K^/2{M\x\) 

(M|x|)3/2 



1/2 



+ 



'm\x\) 



m|x|)^/2 



1/2- 
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As it follows from the above analysis, the full viscosity is given by the formula 

N' 



where the contribution Arj is produced by the interaction between these two types of ex- 
citations. This contribution corresponds to the cross term in the square of the function 
D(x): 



cross term = 24^ ■ 2 



■^^,Ky,{M\x\) Ky,{m\x\) 



1/2 

(20) 



(M|x|)3/2 (m|x|)3/2 
This cross term can be approximated by a function of the type of Eq. (fT3l) as 

cross term ~ 24 ■ N , . , (21) 

{fi\x\yy^ 

where fx = ^^'^"^ . Then, by virtue of Eq. ([H]), Ar] can be expressed through the mass 
parameter fi and the yet unknown coefficient A/" as 

A7] = V8^Af^. (22) 

In terms of the spectral density, this means that the interaction between the two types of 
excitations is also modeled by the Lorentzian Ansatz ([7]), whose width /i is given by the 
mean value of M and m. 

Using now the explicit form of the function -^^^05^, which can be found in Eq. (fTOj) . we 
see that the approximation fl2ip to Eq. (120]) can be written as follows: 



x\ 



To make the right-hand side of this expression really constant, we have to disregard 
in all the three brackets. That is, we restrict ourselves to the leading pre-exponential terms 
in the exponentially falling off functions ^^"^^j/? , ^^[^[^i/? , and ^^0^- This yields the 
following result: 



Owing to Eq. (I22l) . the sought full shear viscosity f lT9l) finally reads 

(23) 



7] 



N N' 2^]^ 
+ ^ + 



M5 m5 ^'^{Mmfl^ 
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At T ^ Tc, the three terms on the right-hand side of Eq. (l23l) . once divided by T^, scale 
(with the double logarithmic accuracy) as 0{g^), 0{l/g'^), and 0{g^), respectively. We see 
that, at such temperatures, Arj is parametrically larger than the contribution of excitations 
with the correlation length 1/M by a factor of 0{l/g^), while being parametrically smaller 
than the contribution of excitations with the correlation length 1/m by a factor of 0{g'^). 

In order to evaluate Eq. (123|) numerically, we adopt the lattice-simulated SU(3) Yang- 
Mills critical temperature = 270 MeV and the two-loop running coupling jol 
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* = 26„ 1„ - + ^ 1„ ^2 lu - j . where = ^. = - . A = 0.104T„ K = 3. 

For the correlation length ^ of the SU(3) Yang-Mills vacuum in the confinement phase, 
entering Eq. (fT6l) . we use the lattice value from Ref. which reads ^ = 894 MeV. Fur- 
thermore, also in the confinement phase, we should choose such a lattice value of the gluon 
condensate G that corresponds to an Ansatz for the two-point correlation function (fT2|) with- 
out the i^-term. According to Ref. [l2|, this value is G = 2.21 GeV^. Next, the entropy 
density, s = s{T), in which units rj is measured, can be obtained by the formula s = dpi^x/dT, 
where we use for the pressure piat the corresponding lattice values from Ref. j6i]. The value 
~ 2.2 is fixed by imposing cancellation of the latent-heat contribution to the pressure. 
Namely, this contribution is expressed by the discontinuity oi e/T^ at T = Tc, which, ac- 
cording to Fig. 7 of Ref. jsl, is ~ 2.2. In Fig. [H we plot the resulting s/T'^ as a function of 
T/T,. 

Implying the subtraction of the vacuum contribution, and noticing that coth^^j = 
1.07 ~ 1, we use in Eq. (123|) a parametrization of Gt in the form Gt = G ■ {f^ — 1), where Jt 
is given by Eq. (ITTIl . We find the value of upon the comparison with the best known lattice 



value of the --ratio 13], which reads -L -, = 0.134. Then, the value of comes out to 

s I II s I T=l. 65 Tc ' 

be just ~ 1.65 Tc, and the resulting ^ turns out to be a monotonically increasing function 
with the minimum value at T = T^: -L ^ — 0.081. We notice that this value is only a few 
per cent larger than the lower bound of ^ predicted by the AdS/CFT-correspondence |3|. 
The resulting ^-ratio is plotted in Fig. HJ Also, using this value of T*, we plot in Fig. [3] 
the ratios of the 1st and the 3rd terms on the right-hand side of Eq. ( l23l) to the leading, 
2nd, term. As mentioned above, the first of these two ratios goes as 0{g^^), becoming at 
T = 4.5 Tc as small as only 0.01, whereas the second ratio goes as 0{g'^). 

Furthermore, we evaluate the extent to which the correcting factor ([9]) at |A;| ~ 1 may 
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1 1.5 2 2.5 3 3.5 4 4.5 

T/Tc 

n 

Figure 1: The (s/r^)-ratio corresponding to the lattice values for the pressure from Ref. J6|] (courtesy 
of F. Karsch), as a function of T/Tc. 

differ from 1 for temperatures T G (Tc, T*). It turns out that, regardless of the argument, 
m, M, or /i, entering this correcting factor, the maximum value it can reach is < 1.42. 
Consequently, at T G (Tc, T*), the resulting ^ may be up to 42% larger. That would lead to 
a discontinuity at T = of the curve depicted in Fig. O since in the dimensionally-reduced 
phase of T > the correcting factor ([9]) is always equal to 1. This uncertainty, emerging 
due to possible deviations of the correcting factor from 1 at T G (Tc,T*), is related to a 
nonuniversal character of the Lorentzian Ansatz ([7j). 



III. CONCLUDING REMARKS 



In this Letter, we have considered shear viscosity within a model of the gluon plasma 
based entirely on the nonperturbative stochastic fields which, however, may develop excita- 
tions with two different correlation lengths, ~ {g'^T)^^ and ~ [g'^T)^^. Excitations of the 
first type possess all the features of the genuine nonperturbative stochastic chromo-magnetic 
fields, which survive the deconfinement phase transition in Yang-Mills theory (qI. Rather, 
excitations of the second type are introduced for the sole purpose to reproduce, within one 
and the same model, the known high-temperature perturbative contribution to the shear 
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Figure 2: The evaluated (ry/s)-ratio as a function of T/T^. 

viscosity. Bringing both types of excitations together, we have studied the contribution of 
nonperturbative stochastic Yang-Mills fields to the shear viscosity, relative to the pertur- 
bative contribution, as well as the result of the perturbative-nonperturbative interference 
(cf. Fig. [3]). Remarkably, the interference contribution parametrically exceeds the purely 
nonperturbative one by a factor of 0{1/ g^). Furthermore, using the best known lattice value 
for the ^-ratio from Ref. 13|, we have evaluated the temperature dependence of this ratio, 
which is depicted in Fig. O The minimum value of the resulting ^-ratio, reached at T = Tc, 
turns out to be only a few per cent larger than the lower bound of predicted for this 
quantity by the AdS/CFT-correspondence. 
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Figure 3: The ratios of the 1st and the 3rd terms on the right-hand side of Eq. (|23|) to the 2nd 
term are represented by the lower and the upper curves, respectively. 
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